Statistics of subgap states in s± superconductors 



O 

(N 

PL, 



C ■ 
O 1 

o 
u ■ 

s : 

a ■ 



i 

c 

o 
o 



> 
m 

m 
o 

O 
O 



X 



A. Glatz and A. E. Koshelev 
Materials Science Division, Argonne National Laboratory, Argonne, Illinois 60439 

(Dated: February 2, 2010) 

There is strong support in favor of an unusual s± superconducting state in the new iron-based 
superconductors, in which the gap parameter has opposite signs in different bands. In this case 
scattering between different bands by impurities has a pair-breaking effect and introduces states 
inside the gap. We studied the statistics of disorder-induced subgap states in s± superconductors 
due to collective effects of impurities. Numerically solving the two-band Bogolyubov equations, 
we explored the behavior of the density of states and localization length. We located the mobility 
edge separating the localized and delocalized states for the 3D case and the crossover between the 
weak and strong localization regimes for the 2D case. We found that the widely used self-consistent 
T-matrix approximation is not very accurate in describing subgap states. 



The recent discovery of high-temperature supercon- 
ductivity in the iron arsenide LaFeAsOi-xF^ followed 
by the discovery of several new classes of superconduct- 
ing materials^ led to a major breakthrough in the field of 
superconductivity. The Fermi surface of these materials 
is composed of several electron and hole sheets located 
near different points of the Brillouin zone^, see Fig. QJt. 
There is theoretical reasoning in favor of an electronic 
origin of superconductivity and a unusual superconduct- 
ing state in which the order parameter has opposite signs 
in different bands (s± state) X This scenario is supported 
by the observation of a resonant magnetic mode in the 
superconducting state by inelastic neutron scattering. 5 

In the s±-state interband scattering due to poten- 
tial impurities (see Fig. [T^) has a pair-breaking ef- 
fect similar to magnetic impurities in conventional 
superconductors^— Such impurities introduce states in- 
side the gap known as Shiba-Rusinov states for the 
magnetic-impurities problem^— This leads to a finite 
density of states (DoS) at zero energy and strongly in- 
fluences the superconducting properties at low temper- 
atures. Most iron pnictides are doped superconductors 
and disorder due to dopant atoms is unavoidable. Ex- 
perimental properties of iron pnictides which may be ex- 
plained by a finite DoS at the Fermi level due to pair- 
breaking disorder include the "residual" linear tempera- 
ture dependence of the specific heat in the superconduct- 
ing stated and a quadratic temperature dependence of 
the London penetration depth at low temperatures^ 

Transport properties at low temperatures, such as 
thermal conductivity 12 and microwave surface resistance, 
are sensitive to localization state of the low-energy quasi- 
particles. A random impurity potential may localize 
quasiparticles within some energy range. For 3D super- 
conductors we typically expect two regimes: (i) At small 
concentrations of impurities states at the Fermi level are 
localized and a mobility edge at a finite energy exists, 
similar to impurity bands in semiconductors; (ii) At suf- 
ficiently high impurity concentrations all states are delo- 
calized. The critical impurity concentration depends on 
the scattering properties of the impurities. In the first 
regime, the localized states near zero energy contribute 
to the low-temperature behavior of the thermodynamic 



but not transport properties. 

A standard analytical approach to describe collec- 
tive effects of impurities in superconductors is the self- 
consistent T-matrix approximation (SCTM). For the 
Born limit, it was elaborated in the famous paper by 
Abrikosov and Gor'kov 6 and was later generalized for 
strong scattering 7 . Most theoretical works on collective 
effects of impurities in superconductors are based on this 
approximation, see, e.g., the review While giving a 
qualitative description of the impurity band, the SCTM 
approximation has serious deficiencies. For small impu- 
rity concentrations it predicts a hard gap in the spectrum. 
In reality, the DoS is finite at all energies, it has an expo- 
nential taili^ due to rare fluctuation configurations of im- 
purities, similar to the Lifshitz tail in the impurity band 
of semiconductors. Another deficiency of the SCTM de- 
scription is that it ignores localization properties of the 
states. Localization of quasiparticles in superconductors 
was studied for dirty d-wave superconductors^ and for 
the mixed state of disordered s-wave superconductors 15 . 
To our knowledge, localization in the impurity band of 
superconductors with pair-breaking impurities was never 
studied. 

Several recent papers address different aspects of the 
impurity-induced subgap states^— and their possible 
influence on properties of iron pnictides 19 . Studies of 
collective impurity effects, however, do not go beyond 
the SCTM approach. Motivated by the importance of 
disorder-induced subgap states for the properties of iron 
pnictides and the absence of an accurate theoretical de- 
scription of these states, we performed a detailed study 
of their statistical properties based on numerical calcu- 
lations. We explore the behavior of the density of states 
for s± superconductors as function of scattering parame- 
ters and impurity concentration and compare the results 
with the SCTM approach. We also explore localization 
properties of states in order to locate the mobility edge 
in the parameter space. 

Our study is based on the two-band Bogolyubov equa- 
tions for the two-component wave functions, $ a = (""), 
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FIG. 1. (color online) (a)The schematic band structure of iron 
pnictides. Solid and dotted arrows illustrate intraband and in- 
terband scattering events. (b)Evolution of the 2D subgap DoS 
with increasing concentration of impurities Hi for isotropic 
scattering. Symbol sizes are proportional to L cn f ([3]). 



Here a is the band index, fi are Pauli matrices in Nambu 
space, i a = £ a (k) — ep « v^ Q (k — kjr a ), A Q are the 
gap parameters. We assume that A2 = — Ai. The last 
term in the equation describes the interaction with im- 
purities. The interband scattering is described by the 
off-diagonal terms in the matrix U l a p. For scattering be- 
tween bands 1 and 2, separated by wave vector Q equal 
to half of the reciprocal-lattice vector, U l 12 contains the 
factor exp(iQR ; ) which only takes values ±1 depending 
on R; . This means that even for identical impurities U[ 2 
has random signs and its average is zero. We neglect in- 
homogeneities of the gap parameters due to impurities. It 
is known that these inhomogeneities are small and do not 
influence the quasiparticle states much. The key param- 
eter of an isolated pair-breaking impurity is the energy of 
a localized stated, Eq/A = eq = \/l — 4r c ff, where we 
introduced the effective interband scattering parameter, 
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with j a p = iri/ a U a p being the reduced scattering ampli- 
tudes and v a = j ff^p- £(£a(k)~ £f) being the normal DoS 
(per spin) for band a. 

We explore the properties of the subgap states for 
the two- and three-dimensional cases. For the numerical 
analysis, we rewrite the equations in a form containing 
only wave functions at the impurity sites, 

$ (R,) = ^.g Q (R - R,0f z7 ^(R,0, (2) 

0,1' 

where the reduced Green's function is defined as 
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and d is the spacial dimensionality. For each impurity re- 
alization inside a box of size R d we find the set of eigenen- 
ergies E\ and the corresponding four-component wave 
functions # ; A = [*^(R ; ), *2 (R-0], see appendix[B] From 
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FIG. 2. (color online) Left: Evolution of the subgap DoS with 
changing scattering rate for isotropic scattering for two im- 
purity concentrations. Symbol sizes are proportional to the 
confinement length. Right: Comparison of the DoS shapes for 
the 2D and 3D models with identical pair-breaking parame- 
ters and SCTM results for two impurity concentrations. 



the set of eigenenergies, we compute the average DoS, 
N S (E) — (J2\ ${E — E\)), where the average is taken 
over many impurity realizations. We normalize N S (E) 
to the total normal DoS for excitations, N n = Av. To 
characterize localization properties, we compute the av- 
erage confinement length for given energy, 



L C nf(E, R) 




= EiR^tl 2 («* = 1,2), and 




(3) 



where (r r a 

J2 a [K(R)| 2 + K(R/)| 2 ]- The behavior of L cnf (E,R) 
with increasing system size, R, determines the nature of 
the states. For delocalized states L cn {(E, R) is limited by 
the system size and linearly grows with R. For localized 
states L cn f(E, R) saturates at a finite value, which gives 
the average localization length of states with energy E, 
hm^^oo L cn{ (E, R) — L\ oc (E). 

We study the subgap densities of states at different 
concentrations of impurities and scattering parameters. 
We consider first the case of isotropic scattering, when all 
matrix elements are equal. Figure [T}) shows the evo- 
lution of the subgap DoS with increasing concentration of 
impurities rii for a moderate scattering rate, = 1 for 
all a and (3, and the system size R = 5£. The energy is 
normalized to the gap value A and the coherence length 
is defined as £ = Vp/A.. For small concentrations of im- 
purities the DoS has a peak near the energy of the local- 
ized state. With increasing impurity concentrations, the 
peak broadens and becomes completely smeared already 
at relatively small impurity concentration n^ 2 =2. At 
higher concentrations the DoS becomes almost constant 
and comparable with the normal DoS. The left part of 
Fig. [5] shows the subgap DoS for two impurity concen- 
trations for a wide range of isotropic scattering rates. 
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FIG. 3. (color online) Plots (a) and (c) show the dependences 
of the zero-energy DoS on the scattering strength for isotropic 
scattering for 2D and 3D cases. Solid lines show the SCTM 
results. Plots (b) and (d) show the size dependences of the 
confinement length at different scattering strength. 



The energy location of the peak for small concentrations 
goes down with increasing scattering strength, while the 
maximum DoS is almost independent on y a /3. The peak 
smears with increasing impurity concentration and the 
magnitude of the density of states at large concentra- 
tions depends only weakly on the scattering strength for 
7a/3 > 0.5. Also, for large scattering rates the DoS ap- 
proaches a limiting shape corresponding to the unitary 
limit. For small scattering strength, — 0.2, we found 
the typical "Lifshitz tail" behavior. 

Within the SCTM approximation (see appendix [A"| the 
DoS does not depend on dimensionality of superconduc- 
tor, it is determined by the location of the single-site en- 
ergy level and the pair-breaking parameter proportional 
to the impurity concentration a = 2niT e g/(irisA). Even 
though the overall evolution of the DoS shape with vari- 
ations of the scattering strength and impurity concentra- 
tion is similar in the 2D and 3D cases, the universality 
suggested by the SCTM approach does not exist. In the 
right part of Fig.[5]we compare the computed DoS for two 
impurity concentrations corresponding to the same pair- 
breaking strength in the 2D and 3D cases with the SCTM 
results. One can see that the DoS shapes for the two 
dimensionalities are similar but not identical. One no- 
ticeable difference is that the 3D DoS is typically smaller 
at low energies. The SCTM approximation does not re- 
produce the DoS shape at low impurity concentrations, 
the sharp peak in the center and the small features at 
the sides are not reproduced. These features appear due 
to the oscillating dependence of the impurity pair energy 
on their separation and correspond to the pairs separated 
by the distance at which the energies have extrema (see 
appendix [C]). For large impurity concentration we found 
a pronounced dip in the DoS for energies slightly smaller 
than A, also not reproduced by the SCTM approxima- 
tion. 



We now focus on the region near zero energy, which 
determines the low-temperature behavior of the super- 
conducting parameters. Figures 03a, c ) show the depen- 
dences of the zero-energy DoS on the isotropic scattering 
strength for fixed impurity concentration for the 2D and 
3D models. We see that the DoS has a negligible size 
effect even at the smallest studied sizes. The SCTM ap- 
proximation only roughly reproduces the the shape of 
the numerical curves and does not describe the tail re- 
gion. The size dependences of the confinement length are 
shown in Fig. EIb,d). We can see that for the 2D case 
at 7q/3 < 0.6 the confinement length saturates at large 
R approaching a finite localization length. No clear sat- 
uration of L cn { is observed for 7 Q( g > 0.6. This may im- 
ply that states are delocalized or the localization length 
may be much larger than the studied system sizes. The 
second possibility looks more plausible, because for 2D 
disordered systems all electronic states are expected to 
be localized. The region rs 0.6 probably marks a 
sharp crossover between the weak and strong localization 
regimes. We actually observe a noticeable downward cur- 
vature in the dependences of £ C nf vs. R for all 7 Q /3 which 
may be interpreted as a tendency towards localization at 
larger length scales. For the 3D case we expect a true 
mobility edge which we can estimate as the value of 'jap 
at which L cn f has a clear saturation tendency at large R. 
The estimated location of the mobility edge, 7 Q( g ps 1.25, 
is marked in Fig. [3p. It is close to the critical value 
bounding the SCTM gapped regions. We found that the 
DoS values at the mobility edge are quite small, (~ 0.02 
in our example). These values are significantly smaller 
than the DoS at the localization crossover for the 2D case. 
This observation implies that for identical bands in the 
3D case there is only a very narrow parameter window 
within which the localized states at zero energy provide 
a noticeable DoS. 

In figure [¥Ja,b) we compare the numerical and SCTM 
dependences of the zero-energy 2D DoS on the impuri- 
ties concentration for weak and strong isotropic scatter- 
ing. We also show the concentration dependences of the 
transition temperature evaluated using the Abrikosov- 
Gor'kov formula (see Ref. [f| and appendix For small 
scattering strength a noticeable DoS appears at the Fermi 
level only when T c is strongly suppressed. In contrast, for 
large scattering strength the DoS reaches values compa- 
rable with the normal DoS already at very minor sup- 
pression of T c . 

In Fig. 2J: we present the dependences of the zero- 
energy DoS on the scattering strength for three concen- 
trations of impurities, n^ 3 = 1, 2, and 4 and compare 
them with the predictions of the SCTM approximation, 
which only coarsely reproduces the evolution of the DoS 
with increasing concentration. It does not describe the 
tail regions and systematically underestimates the value 
of the zero-energy DoS in the unitary limit, correspond- 
ing to large j a p. 

Up to this point we considered the case of isotropic 
scattering when all scattering amplitudes are equal. In 
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FIG. 4. (color online) Plots (a) and (b) show the dependences 
of the zero-energy 2D DoS and transition temperature on con- 
centration of impurities for two isotropic scattering rates. The 
clean-limit coherence length is taken as /cf£ = 2.63. (c)The 
dependences of the zero-energy 3D DoS on the scattering 
strength for three concentrations of impurities. (d)Densities 
of states as function of the interband scattering parameter 
r e ff, Eq. ([T]), for two simulations series, isotropic scattering, 
711 = 712, and fixed intraband scattering amplitude, 711 = 1. 



general, one can expect that the interband scattering 
is always weaker than the intraband one. Within the 
SCTM approach the DoS is not sensitive to the individ- 
ual scattering amplitudes and is completely determined 
by the parameter r e fj, Eq. ([1]). To clarify the role of 
the relative strength of the scattering amplitudes, we 
study the subgap states for different interband ampli- 
tudes 712 = 721 at fixed intraband amplitudes 711 = 722- 
The DoS dependences on T e g for the two simulation se- 
ries are compared in Fig. |4ji. We can see that the DoS is 
not simply determined by r c ff, as suggested by the SCTM 
approach, but also sensitive to the relative strength of the 
intraband and interband scattering. For the same T c ff the 
DoS decreases with decreasing ratio 712/711- 

In conclusion, we explored the subgap DoS and lo- 
calization properties for disordered s± superconductors. 
We found that the widely-used analytical description 
(SCTM) is incomplete and not very accurate. Disorder 
makes superconductivity "gapless" , the DoS at E = is 
always finite. In the 3D case there is a mobility edge sep- 
arating localized and delocalized states. It reaches zero 
energy at a critical impurity concentration above which 
all states become delocalized. In the 2D case the mobility 
edge is replaced by a crossover separating strongly and 
weakly localized states. The development of quantitative 
theory of the subgap states is crucial for the understand- 
ing properties of the iron pnictides and other supercon- 
ductors with pair-breaking impurities. 
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Appendix A: Self-Consistent T-matrix 
approximation 

The self-consistent T-matrix approximation^ is de- 
fined by the coupled equations for the average Green's 
function G a p(k) = SapGaCk) and the self-energy part 



G„(k) 



E - i a f z + A a r x 



(Ala) 



l-T z U a pG o , ) r z U 0a ) (Alb) 



with Go, a = J (2^ (k) . Note that the off-diagonal 

components of Y> a /3 average to zero due to the random 
signs of the off-diagonal components of U a p . The partial 
density of states is related to the Green's function as 

N a (E) = Im |Tr Go, a (E — iS) j. Using the expansion 

S a = J2i ^Q.i^i with i = 0, x, y, z, we can express N a (E) 
via the ratio u a = (E — S Q ,o) / (A — S aj2: ) as 



N a (E) = 2v a \m 



(A2) 



Here the factor 2 accounts for the electron and hole ex- 
citations. 

In the case of identical bands, this system reduces 
to an equation for only one complex parameter u = 
(E — S ) / (A — S x ), similar to the Shiba equation for 
magnetic impurities^, 



u 1 



a- 



E 
A : 



_ 2n a r Gff 

TTl/A 

2n, 



n 2 A 



(A3) 
(A4) 
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7wA 1 + 7| 2 + 7^ + 2712721 + (722711 - 712721) 2 

We observe that increasing the intraband scattering po- 
tential 7 QCe increases the intraband scattering time T12 
and diminishes the pair-breaking parameter a,*^ We also 
note that the SCTM results do not depend on dimension- 
ality of the superconductor. 

We summarize the most important analytical results 
obtained within this approximation. The SCTM approx- 
imation gives a gapped state for a < e\ and gapless state 
for a > Sq with finite total density of state at zero energy 



A^(0) = Au 



V" 2 - 4 



(A5) 
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We stress again that the existence of a gapped state is 
unrealistic feature and deficiency of this approximation. 

The suppression of T c is determined by the Abrikosov- 
Gor'kov formula^, 



Here An is the gap parameter for the clean case. 



Appendix B: Numerical simulations 



In |^=^(l/2 



l/27rn 2 T c ) - V(l/2), 



(A6) 



where ip( x ) ls the digamma function. This famous result 
is almost always understood too literally. In fact, it just 
gives an approximate typical value of the critical temper- 
ature. In real systems, due to the random arrangement of 
impurities, the transition temperature is inhomogeneous 
and the transition has percolative nature. 

The average zero-temperature gap parameter is deter- 
mined by the equations^ 



hr^ = -- 
A 2 1 



for a < e n 



so 



(A7a) 



To develop a precise theoretical description of the sub- 
gap region, we solve Eqs. (2) of the main paper nu- 
merically for the two-dimensional and three-dimensional 
cases. The main element of these equations is the Green's 
function in real space, g a {r)- At r — the Green's func- 
tion does not depend on dimensionality, 

G 1 i A S- 

9c(0) 



Large-r asymptotics of g a (r) for kp, a r 3> 1 are given by 



with a = uqAq/A and 



2e + l +h \ UQ+ V 1 + U °)-^+4 



p 2 -l 



with Uq 



, U 

arctan uo — £o arctan — 



(A7b) 



a 



-0' 



for a > e 



Ja(r) 



-E + A a t x 



(kF, a r - - t z sin (k F , a r - 

exp(-y/Al-E 2 r/v F , a ) 



sin(fci?. Q r) — f z cos(kF, a f) 



kF, a r 



cxp (-y/Al- E*r/v F , a 
for 3D case 



for 2D case 



As the probability to find two impurities at distance 
~ l/kp is very small, the structure of states is mostly 
determined by these asymptotics and the value at r = 0. 
To match the large-r asymptotics with the r = value, 
we use an approximate forms of the Green's functions. 
For the 2D case we use 



-E + A a r x 



Jo (kF, a r) + t z Ji (k F , a r) 



expl-y/Al-EPr/vF,*) , 



(Bl) 



where Jo (x) and J\ (x) are Bessel functions and for the 
3D case we use 



E+A a f x sui(k F , a r) A .... cos(k F . a r) 
- T z U{k F . a r) 



y/Al-E 2 k F , a r 
exp(-\/A| - E 2 r/v F . 



(B2) 



with the interpolation function U{z) = z 2 /(z 2 + 1). We 
expect that the exact behavior of the Green's functions 
at A;ir jQ ,r~l not reproduced by these interpolations have 
little influence on the properties of the subgap states. In 
this paper we limit ourselves with the simplest case of 
two equivalent bands, meaning that kp.i — kp.2 = kp, 
Vf,i = v f.2 = vf, and Ai = — A 2 = A. The gap A 
is used as a unit of energy and kp 1 is used as a unit of 
length. 
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We employ the following numerical procedure. First, 
we define an impurity realization by Ni random coordi- 
nates [Hi] in a box, < Ri. a < R, and random impurity 
signs Si = ±1 for the off-diagonal scattering amplitudes, 
U[ 2 = 5iU\2- From the linear 4iVj x 4iV; system defined 
by Eq. (2) of the main paper, we find the eigenenergies, 
E\, and corresponding eigenstates V I' A (R;). From the set 
of eigenenergies, E\, we compute the average density of 
state 



2D, 7a = 1,^=5 
n.^ 2 =0.4 




N S (E) = - E x ) 



(B3) 



where average is taken over many impurity realizations. 
Practically, this implies that for every realization we 
find the number of states AN(E) falling within the en- 
ergy interval [E — AE/2, E + AE/2] and then compute 
the average N S (E) = (AN{E)/(AER 2 )}. As an iso- 
lated impurity generates one localized state, for small 
concentration of impurities the normalization condition 
J A N s (E)dE = Ni/R 2 is satisfied. We normalize N S (E) 
to the total normal density of states for excitations, 
N n = 4z/, where for the 2D case the single-band DoS 
per electron is given by v — kp / {2nvp) and for the 3D 
case, v — k 2 F / '(2ir 2 vp). 

To characterize localization properties, we also com- 
pute the average confinement length for states at given 
energy, 



L cn{ (E,R) = { M 



(B4) 



E X =E 



where a — x,y[, z] and the confinement length of state A, 
(^ r a) v ^ s determined by its wave function VE^ as 



A ' 



' •• A = y^,Rl,a\^l ' 2 



i 2 , (4 = EC 



A I 2 



,A|2 



with |^| 9 = E Q [l^(R/)| 2 + l^(R/)| 2 ] and EJ^, , 
1. The behavior of the confinement length with in- 
creasing system size, R, determines wether states at 
given energy are localized or not. For delocalized states, 
L C nf(E, R) is limited by the system size and grows pro- 
portionally to R. For localized states, L cn f(E, R) satu- 
rates at a finite value, which gives the average localization 
length of states with energy E, lirnft_ J . 00 L cn { (E, R) = 
L\oc(E). In three dimensions for relatively small con- 
centration of impurities one can expect the existence 
of a mobility edge in the subgap region separating lo- 
calized and delocalized states. It is defined as the en- 
ergy EME^ni, Jap) a t which the localization length di- 
verges, L\ oc {E — > Eme) ~ > oo. At a critical concentra- 
tion of impurities depending on the scattering param- 
eters, n cr {^j a p), the mobility edge reaches zero energy, 
EuE[n C r {la.fi) j Ja.p] = 0, and all states become delocal- 
ized. 



FIG. 5. (color online) Left: The 2D DoS for small concentra- 
tion of impurities, same as in Fig.l of the main paper. Right: 
The energy of impurity pair as a function of the distance be- 
tween impurities for same-sign and opposite-sign impurities. 
The arrows mark the extremum points which account for the 
small peaks in the DoS. 



In the calculations of the concentration dependences 
presented in the Figs. 4a, b we took into account sup- 
pression of the average gap parameter described by Eqs. 
(|A7a[) and (|A7b[) and corresponding increase of the co- 
herence length £ = vp/A. 



Appendix C: Density of states at small 
concentration of impurities: pairs-dominated regime 

At small concentrations of impurities, rij£ d <§; 1, the 
density of state is determined by impurity paired The 
interaction between two close impurities gives the correc- 
tion to the energy^, E p ^ T (r) = So + 5e(r) and the behav- 
ior of the correction 5e(r) depends on the relative sign 
of the off-diagonal scattering potential XJ\i for two im- 
purities. For same-sign impurities the single-site energy 
level splits into two levels corresponding to symmetric 
and antisymmetric combinations of the wave functions 
at the impurity sites. The separation-dependent energy 
corrections, Se±, rapidly oscillate with distance between 
the impurities r as 



sm(k F r + a d ) 
fe±(r)0C± (fc f r)C-i)/» 6XP 



where £ = Vp / A is the coherence length and d is space di- 
mensionality, see Fig. [SJ For opposite-sign impurities the 
energy level remains double-degenerate with two states 
corresponding to localization near two impurity sites. 
The separation-dependent energy shift in this case is al- 
ways positive and much smaller than for the same-sign 
impurities, 



Se(r 



sin (kpr + ad) 
oc — - — — j ; — exp 



(kpR) 



d-l 



-2J1 -efrlt 



The coefficients in the energy corrections depend on the 
scattering parameters of impurities. 

The contribution to the DoS coming from impurity 
pairs with separation r p less than the typical distance, 
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r p <C n~ , can be evaluated in a simple way^ The con- 
centration of the impurity pairs with separations between 
r p and r p +dr p is given by ^r? i r^ 1 ~ x dr p with A 2 = 2ir and 
A3 = 47r. The contribution to the DoS at the energy E is 
given by the pairs satisfying the equation E = E va \ r (r p ) 
and for a nonmonotonic dependence -Epair (**), this equa- 
tion may have several solutions. The pair contribution 
to the DoS can be evaluated as 

N S (E) = ^n^[r p (E)]^ ^ , 

where the sum is taken over all values of r p corresponding 
to the same energy. 

For pair separations corresponding to the energy ex- 



trema, E pa it(r) « E e + §(r — r e ) 2 , the isolated-pairs ap- 
proximation gives a divergency at E — > E ei 

Ns(E) „ Adk£L. 

^/2a (E - E e ) 

These singularities, smeared by interactions with more 
remote impurities, account for the small peaks found in 
the DoS at small concentrations, see Fig. [5] The en- 
ergy for opposite-sign impurities has a series of minima 
at exactly E = Eq meaning that there is also such a pair 
singularity at the peak center. This explains the sharp- 
ness of the peak at small concentrations of impurities. 
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